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We derive an analytic formula which describes the final bound state dependence in ultracold 
three-body recombination. Using an energy-dependent loss parameter, the recently observed broad 
resonance in an ultracold gas of 6 Li atoms [1, 2] is described quantitatively. We also provide an 
analytic and approximation for the three-body recombination rate which encapsulates the underlying 
physics of the universal three-body recombination process. 
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Over the last several years, the field of Efimov physics 
has generated a great deal of excitement [3, 4]. Ultracold 
atomic gases provide a unique environment for studying 
these exotic states. Through the use of a broad Fano- 
Feshbach resonance the s-wave scattering length can be 
tuned over several orders of magnitude [5-7] . Using these 
tools many universal behaviors associated with Efimov 
physics have been observed in both three- and four-body 
recombination processes [8-12], largely confirming many 
theoretical predictions (see Refs. [13-16] for some exam- 
ples). The majority of these experimental studies have 
taken place using an ultracold gas of bosonic atoms. 

Recently, a new class of Efimov state has emerged in 
a three-component degenerate Fermi gas of 6 Li atoms. 
For magnetic fields between 30 and 600 G, two reso- 
nances in three-body recombination have been observed 
[1, 2]. The positions of these resonances corresponds to 
a single three-body bound state crossing the continuum 
at lower and then again at higher field strength. While 
the position of these resonances is fit very well by exist- 
ing three-body recombination theory [17-19], the widths 
are not. The standard theory, in which it is assumed 
that the short range behavior of the three-body system is 
largely independent of the strength of the external mag- 
netic field, predicts that the two resonances should have 
roughly the same width. The experimental evidence show 
a higher field resonance that is significantly broader than 
the first. 

Three-body recombination occurs when two particles 
combine to form a dimer state releasing the resulting 
binding energy in the form of kinetic energy between 
the dimer and a third particle. It has been proposed 
by Wenz et al. [20] that the binding energy Ef, of the 
final dimer state is to blame for the observed discrep- 
ancy in resonance widths. The binding energy of the 
first several deeply bound two-body states in 6 Li have a 
strong magnetic field dependence [20]. By assuming a 
1/Ef, dependence in the loss parameter they were able 
to find reasonable agreement with both resonances. The 
question still remains, however, as to where this energy 
dependence comes from and whether l/E^ is the correct 
form. 

In this paper we present a simple mechanism that de- 



scribes the broadening of a three-body recombination res- 
onance due to the magnetic field dependence of a deeply 
bound two-body state. We proceed by assuming that this 
precess can be described by a two-channel inelastic scat- 
tering process within the adiabatic hyperspherical repre- 
sentation, shown schematically in Fig. (1). The upper 
channel consists of three free particles scattering near 
threshold, while the lower channel consists of a tightly 
bound dimer with binding energy A and a free particle. 
The two channels are coupled at some small hyperradius 
ro by a non-adiabatic coupling P-matrix element: 



P if (R) = <*/ (i?; n) \d/dR\ $>* (R; Q)) s 



(1) 



where 3>j(/) (R~, O) is the hyperangular channel function 
describing the initial (final) adiabatic channel and the 
matrix element is taken over the hyperangular degrees of 
freedom. 




FIG. 1: (color online) A Schematic of the two-channel scatter- 
ing process that controls threshold three-body recombination 
is shown. 

Mchta et al. [16] showed that the inelastic scattering 
rate at the three-body threshold for this type of system 
can be described within the Wentzel-Kramers-Brillouin 
(WKB) approximation as 



9h 5 Ce~ 2 ' 1WKB sinh 2-q 
4m 3 E 2 cos 2 <I>wkb + sinh 2 r\ ' 



(2) 
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where m is the atom mass, E is the initial asymptotic 
kinetic energy of the three atoms, and the constant C is 
determined below. In Eq. (2) jwkb is the WKB tunnel- 
ing integral and 4>wkb is the WKB phase accumulated 
in the inner potential well of the upper potential, i.e. 



Iwkb = Re 



>WKB 



= Im 



pR2 



2 £( Vi (R)-E) + -^ dn 



(3) 



(4) 



where Vi (R) is the initial adiabatic potential which can 
be found using known methods (see for instance Ref. 
[21]) and i?i is the outer classical turning point, and E is 
the initial asymptotic kinetic energy of the three atoms. 
The extra 1/4R 2 in Eqs. (3) and (4) is due to the Langer 
correction. The proportionality constant C ~ 300 in Eq. 
(2) is extracted by comparing the K 3 in case when all 
three scattering lengths are equal to a, where \a\ 3> r , 
to the known recombination rate for this case from Ref. 
[22]: 



K* = 



M590 



a 4 sinh 1r\ 



m cos 2 4>wkb + sinh rj 



(5) 



In Eq. (2), 77 is an imaginary phase which parametrizes 
the losses from the initial channel. The ratio of the out- 
going to incoming probabilities in the initial channel [22] 
is given in terms of r\ as 



out I 



I*,; 



= e' 4 ". 



(6) 



The constants 77 and r are determined by short range 
physics and are usually used to fit the position and width 
of an Efimov resonance [17-19, 23]. For three-body re- 
combination of distinguishable 6 Li atoms, assuming that 
77 and tq are fixed over the experimental range of mag- 
netic fields, there is excellent agreement with the first 
three-body resonance [17-19] at approximately 130 G. 
A problem occurs however at the second resonance; the 
theory prediction is far narrower than the experimental 
results. 

By assuming that the P-matrix element coupling the 
two channels is Lorentzian which peaks at R = r , the 
energy dependence of 77 can be extracted. By consider- 
ing the Landau-Zener transition probability, Clark [24] 
has shown that in this case the probability of making 
a non-adiabatic transition as the system passes through 
the transition region is given by 



1 A 

fiv 8P m£ 



(7) 



where v is the characteristic velocity of the system at 
R = ro, A is the energy separating the initial and final 
channels and P max is the maximum of the P-matrix cl- 
ement. The velocity v is determined by the short-range 
physics of the incoming channel while P max will be de- 
termined by the small R behavior of both the incoming 
and outgoing channel function $i and <&/. 

While P max cannot be determined exactly, we can 
make certain statements about its behavior. We assume 
that the short-range nature of the incoming channel is 
independent of the magnetic field over the range we con- 
sider. The outgoing final channel will be shifted overall, 
but the size and wave function of the final deeply-bound 
dimer state is relatively unaffected. Thus, both P max and 
v are independent of the final state energy A. 

By assuming that any non-adiabatic transition results 
in a three-body recombination event and the system must 
pass twice through the transition, on the way in and on 
the way out, the probability of remaining in the initial 
channel can be extracted: 



I*. 

1 ?. 



= (1 - Pn 



(8) 



Comparing Eqs. (6) and (8) gives a simple equation that 
can be solved for 77: 



1 - exp (-/3A) 
1 



(9) 



hv 4P„ 



With this formula, the unknown short-range dependence 
of r\ is encompassed in a single parameter (3 which is in- 
dependent of the binding energy A. It is important to 
note that we have not increased the number of fitting 
parameters, we have merely shifted the fitting to the sys- 
tem dependent parameter j3 instead of 77. In cases where 
the final state energy does not strongly depend on the 
magnetic field, this extra dependence is not needed and 
77 can be used directly. However in cases where the fi- 
nal state has a strong field dependence, such as 6 Li, the 
proposed parametrization of 77 is more appropriate. In 
general a smaller binding energy leads to a larger loss 
parameter in Eq. (9), and a larger loss parameter leads 
to a broader resonance. Qualitatively, while the binding 
energy dependence in Eq. (6) is similar to the 1/A de- 
pendence assumed in Ref. [20], the detailed behavior is 
considerably different. 

Figure (3) shows the recombination rates predicted 
from Eq. (2) in the threshold regime (incident energy 
E = 10~ 12 a. u.) compared to the experimental results 
of Ottenstein et al. [1] (red circles) and Huckans et al. [2] 
(blue squares). The dashed curve is found by assuming 
that r/ is completely independent of the magnetic field 
strength. The solid curve is found by using Eq. (9) with 
the lowest dimer binding energy from Fig. (2b). Any 



3 




B (in gauss) 



FIG. 2: (color online) (a) The three s-wave scattering lengths 
for 6 Li are shown in atomic units as a function of magnetic 
field, (b) The four smallest dimer binding energies from Ref. 
[20] are shown in atomic units as a function of magnetic 
field. The total nuclear angular momentum and projection 
are given, but are unimportant for the results in this work. 



FIG. 3: (color online) The recombination rate constant pre- 
dicted from Eq. (2) is shown for fixed rj (dashed curve) and 
for energy dependent rj (solid curve) as a function of mag- 
netic field compared to experimental data from Ref. [1] (red 
circles) and Ref. [2] (blue squares). 



one of these four binding energies could have been used 
(with appropriate modifications to /3), but because they 
all have similar magnetic field dependence, the results are 
nearly identical. The fits were found by setting tq = 22 
bohr such that the the first resonance occurs at 130 G. 
The dashed curve was found by setting rj = 0.05 to fit the 
width of the first resonance at B = 130 G using the data 
from Ref. [2]. The solid curve was found by choosing j3 
so that -q — exp(-fiA) — 0.05 at the same B. Both of 
the predictions in Fig. (3) do an excellent job of describ- 
ing the first resonance, while using the A dependent loss 
parameter from Eq. (9) is in astonishingly good agree- 
ment with Huckans et al. [2] and fairly good agreement 
with Ottenstein et al. [1]. 

The initial hyperradial potential used here was derived 
assuming that all scattering lengths were much larger 
than any short range parameters. The hyperradial po- 
tential used in Eqs. (3) and (4) were found by assuming 
zero range interactions, which is appropriate for scatter- 
ing lengths much greater than the size of the two-bod 
interaction. In the case of 6 Li, this size is given by ap- 
proximately the Van der Waals length r$ w 30 bohr. 
Because the smallest scattering length is not too much 
larger than this, we might expect small corrections to this 
due to non-universal behavior in the potentials. Even 
with this caveat, the qualitative agreement seen in Fig. 
(3) is remarkable. 

While the treatment above does give good agreement 
with experimental data, it not necessarily convenient for 
quick comparison to experiment. Examining the three 
scattering lengths from Fig. (2a) shows that they dif- 
fer by a factor of roughly 2 throughout the region in 
which the experimental data is taken. With this, we as- 



sume that they are different enough to use the results 
of Ref. [23] which gives the incoming three-body hyper- 
radial potential for three distinguishable particles with 
large s-wave scattering lengths in four regions: 



Vi (R) = { 
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r < R < h 



(10) 



Here aj, a m and a s are respectively the largest, second 
largest and smallest scattering lengths in magnitude, and 
.so = 1.006 and s% = 0.414 are parameters which are de- 
termined by the universal potential in the limits where 



R < h 



h| and |a s | < R < \a m , h respec- 



tively. It is not clear that these regimes exist for the 
scattering lengths shown in Fig. (2a) , but we will proceed 
assuming that they do. By assuming that transitions be- 
tween different universal behaviors have no significant 
contribution, the tunneling suppression and phase accu- 
mulation can be found using Eqs. (3) and (4): 

<Pwkb = s ln(|a m | / |a s |) + si In (|a s | /r ) , (12) 
s = 1.006; si = 0.414, 

Inserting this into Eq. (2) yields an analytic formula for 
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FIG. 4: The numerically calculated recombination rate con- 
stant from Fig. (3) (solid curve) compared to the analytic 
approximation from Eq. (13) is shown as a function of mag- 
netic field strength. 



energy produces broader resonances when other param- 
eters are held fixed. We have also given an analytic ex- 
pression for the three-body recombination rate in a three- 
component degenerate Fermi gas of 6 Li atoms which pro- 
vides reasonable agreement with the more complex nu- 
merical results. The energy dependent loss parameter 
was derived here using the intuition gained from the adi- 
abatic hyperspherical method, but it is not limited to this 
approach and could be used wherever N-body recombina- 
tion processes are parametrized by an imaginary phase. 
The process of determining the width of loss resonances 
based on knowledge of the dimer binding energy might 
also be turned around. One could use experimentally de- 
termined resonance behavior to find the field dependence 
of deeply bound final states. 
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the three-body recombination rate: 
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Figure (4) shows the results of this formula compared 
to the numerically obtained results in from Fig. (3). 
Both rates use the same A-dependent loss parameter, 
while we have set tq = 32.6 bohr is Eq. (13) to fit 
the resonance position at B = 130 G. The constant 
C univ = 1379 has been c hosen so that Eq. (13) agrees 
with the numerical results at the peak of the resonance. 
While the agreement may not be perfect it indicating 
that the analytic expression from Eq. (13) encompasses 
the majority of the relevant physics. 

In summary, we have derived a simple expression that 
gives the loss parameter 77 as a function of the binding 
energy of the final outgoing states. While this work fo- 
cused on the three-body recombination resonances found 
in a three-component degenerate Fermi gas of 6 Li atoms, 
the method can be applied to other three-body systems. 
For instance in the case of three-body recombination to 
weekly bound dimers, the P-matrix and characteristic 
velocity v are known, and Eq. (9) could be used to 
give a more complete description of the process. This 
parametrization of r\ introduces no extra fitting parame- 
ters and gives excellent agreement with experiment. We 
have shown how a universal Efimov resonance can be 
broadened by the a field dependent binding energy of the 
outgoing final state. In general, a smaller dimer binding 
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